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Abstract 

Using radiative Z° — > t + t~'j events collected with the OPAL detector at LEP at 
^Js = Mi during 1990-95, a direct study of the electromagnetic current at the 
T7 vertex has been performed in terms of the anomalous magnetic form factor 
F 2 of the r lepton. The analysis is based on a data sample of 1429 e + e~ — > 
T + r~'-f events which are examined for a deviation from the expectation with 
F% = 0. From the non-observation of anomalous r + r _ 7 production a limit of 

-0.068 <F 2 < 0.065 



is obtained. This can also be interpreted as a limit on the electric dipole form 
factor F 3 as 

-3.8 x 10~ 16 ecm < eF 3 < 3.6 x 10~ 16 ecm. 
The above ranges are valid at the 95% confidence level. 
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Introduction 



Measurements of the anomalous magnetic moment of the electron ]lj and the 
muon by spin precession methods are considered the most precise tests of 
Quantum Electrodynamics (QED) and are usually expressed in terms of a devi- 
ation of their respective (/-factors from the value of two [Q 
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1159.652193 ± 0.000010) x KT 6 , (1) 



= -) = (1165.9230 ±0.0084) x 10~ 6 . (2) 



Due to the r lepton's short lifetime of (291. Oil. 5) x 10~ 15 s, its anomalous 
magnetic moment cannot in practice be measured by a spin precession method 
and no direct measurement of a T exists so far. While the hadronic and weak 
contributions to a e are very small, they are no longer negligible for a M and a T . A 
theoretical prediction for a T , based purely on QED, is (1173.19 ± 0.01) x 10~ 6 |4[]. 
Additional weak and strong contributions j|, |5| modify this to (1177.3 ± 0.3) x 
10~ 6 . Using the total width of Z° — > t + t~, ref. || indirectly derives an upper 
limit on a T of |a T | < 0.01 at 95% confidence level. 




(P) 



y(k) 



In order to constrain a T as suggested by Grifolz and Mendez 0, we have studied 
the process e + e~ — > r + r _ 7 in which a final-state photon is radiated from one of 
the tau leptons, as shown in the Feynman diagram above. The electromagnetic 
current of a fermion with mass m and charge e can be written using the general 
form factor decomposition 



Jem 



eu{p) 



YF 1 (q 2 ) + —F 2 (q 2 )a^q u + ^ q v F,{q 2 ) 



u(p), 



(3) 



with p',p being the four-momenta of the r lepton before and after the emission 
of the photon with four-momentum q and q 2 = (p — p') 2 . At q 2 = 0, -Fi(O) = 1 
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while F 2 (0) = a/, and eF 3 (0) = df define the anomalous magnetic and electric 
dipole moment, respectively. Note that the Standard Model predicts F 3 = 0. 

However, this ansatz is not directly applicable to the r' — > vertex in 
e + e~ — > r + r _ 7 , since the r', which emits the photon, is off-shell. Instead, the 
pertinent part of the amplitude must be written as 



i(p"+ m) 

—2 2 le 



^F^p' 2 , q 2 ) + ^F 2 (p' 2 , q 2 ) a^q u + ^a^q u F,(p' 2 , q 2 ) 
2m 



e(n)u(p) . 

(4) 



The photon belongs to the final state, so it is real and therefore F 2 (p' 2 ,q 2 ) is 
measured at q 2 = 0, but averaged over a range of p' 2 from m 2 to (Mz — m T ) 2 . In 
this analysis the minimum value of p' 2 after the event selection is (13 GeV) 2 . 

In this paper we search for an excess in the production of e + e~ — > t + t~^( 
events due to a non- vanishing form factor F 2 (p' 2 , 0) as defined by eq. @, assuming 
that F3 = 0. Differential photon rates are compared to Monte Carlo predictions 
for the standard F\ and the anomalous F2 term. The extracted bound on the 
number of excess events from the F2 term is used to determine an upper limit 
on F 2 averaged over p' 2 . Henceforth this interpretation of F 2 is always implied. 
Conversely, assuming F 2 = a limit on F 3 is obtained. Because the sensitivity 
of this analysis is not sufficient to measure a value of F 2 as small as predicted [|J 
by the Standard Model (SM), the reported results mainly address new physics 
phenomena beyond the SM. Such phenomena may occur in the context of com- 
posite t leptons 0, leptoquark models 0, or in models in which the electroweak 
symmetry breaking is driven by the third quark and lepton generation such as 
top-condensation or top-colour models [T(| . 



It should be noted that the ansatz of eq. (§) can parametrize modifications of 
only the r' — > T7 vertex. Radiative corrections involving both final-state taus, as 
well as the non- vanishing p' 2 , therefore limit a direct interpretation of F 2 in terms 
of the r-lepton's anomalous magnetic moment a T = F 2 (0, 0). For physics beyond 
the SM at an interaction scale A new ^> Mz, however, there is no such limitation 
in the above ansatz. In fact, as long as \p' 2 — m 2 \ -C A 2 ew , equating F 2 with the 
a T pertaining to the new interaction is correct. 

The calculation which is used here to predict the distribution of photons 
arising from the different contributions assumes no interference between the F\ 
and the F 2 term. The interference term is suppressed by m 2 /M|. No severe 
restriction is imposed by this assumption for the precision of the F 2 measurement 
described below. Modifications of the results due to the interference term are 
treated at the end of the paper. 
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Monte Carlo simulation 



The Monte Carlo simulation of the process e + e~ — ► t + t~7 with F 2 = F 3 = is 
provided by the program KORALZ |l]J including initial (ISR) and final (FSR) 
state photon radiation up to 0(a 2 ). To the extent that the expectation for F 2 
within the SM is small compared to the sensitivity of this analysis, KORALZ is 
assumed to represent the SM expectation throughout this paper. The r decay is 



simulated by the TAUOLA package which includes photon radiation from the 
leptonic decay products up to 0(a) and also from hadronic decay products using 
the program package PHOTOS fI5|| . According to studies using the KORALZ 
MC, the only source of photons contributing to the selected events studied in 
this analysis will be from ISR and FSR. Photons from tt° decays do not enter as 
background to this analysis after the event selection. 

The contribution of t + t~7 events coming from a non- vanishing form factor 
F 2 is simulated using a calculation by Zeppenfeld |15|] based on the F 2 term in 



eq. (|4]), assuming m T = and neglecting interference. The resulting differential 
cross section is given in the Appendix. In fact, the approximation of m T = 
implies a chirality (=helicity) flip in the amplitude for the F 2 contribution, while 
the Standard Model radiation always conserves chirality. As a result, there is no 
interference between the Standard Model and the F 2 contribution in the massless 
limit. Conversely, the size of the interference term then checks the validity of 
the m T = approximation. A very recent calculation (|T(J of radiative tau pair 



production through anomalous electromagnetic couplings including interference 
effects and a finite r mass confirms the validity of the assumptions (m T = 0, 
interference neglected) made here. Ref. |16[ concludes that anomalous contribu- 



tions from initial-state final-state interference, Z°/7 interference and 7 exchange 
can also be safely neglected. 

Events generated from both the Fi bremsstrahlung term (KORALZ) and the 
F 2 contribution are processed through a full simulation of the OPAL detector 
[ |I1}] . For the purpose of the efficiency determination for the F 2 contribution 
(signal), events have been generated by KORALZ and selected according to the 
5-dimensional differential F 2 cross section (see Appendix) employing a 'hit or 
miss' method. 



Figure H] shows comparisons of the anomalous contibution and the KORALZ 
prediction in simulated distributions of the energy of the radiated photon E 1 (a), 
the acollinearity angle 6^01 between the r lepton directions (b), and the emission 
angle of the photon with respect to the beam direction cos Q 1 (c). Note, that 
the anomalous part is arbitrarily normalized. The striking difference between 
the distributions suggests that these variables are useful discriminators for this 
analysis: F 2 photons appear to be preferentially at high energies and are emitted 
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Figure 1: Comparison between the F 2 signal and the Standard Model expectation 
represented by KORALZ. (a) the photon energy E 1 , (b) the acollinearity angle 
$acoi of the t leptons and (c) the photon angle cos # 7 w.r.t. the beam direction. 
The relative normalization of the distributions is arbitrary, no detector effects are 
included. The dashed lines in (a) and (c) indicate the acceptance cuts, (d) photon 



energy E~ t 
simulation. 



vs. 



acollinearity angle for the F 2 signal prediction after full detector 
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at large angles to both r's. However, E 1 and are strongly correlated (fig. |T] 
(d)), and the cos 7 distribution is almost isotropic in an accepted 6* 7 range with 
reduced background (| cos 9 7 \ < 0.78). Consequently, the benefit obtained when 
using two-dimensional information in the above variables has been found to be 
marginal and also more sensitive to systematic effects. Therefore, in what follows 
only the photon energy distribution is used. 



For the simulation of background processes the Monte Carlo generators |L7 
JETSET 7.4 (qq), RADBAB 2.0 (e+e"), KORALZ 4.0 (ji+fT) and VERMA- 
SEREN 1.01 (2 7 ) have been used. 



Event Selection 



For this analysis events recorded with the OPAL detector during the years 1990 
to 1995 at a centre of mass energy = Mz, corresponding to an integrated 
luminosity of about 180 pb -1 , have been used. 'Off-peak' data were not used to 
avoid deficiencies due to the lack of Z°/7 interference effects in the signal MC. 
The number of produced r pairs is about 230 000. The OPAL detector and its 



performance have been described elsewhere ||18|| . Isolated final-state photons are 
detected in the lead glass electomagnetic calorimeter covering an angular range 
in the barrel region of | cos 9\ < 0.81 with an energy resolution of &e/E ~ 



12%/yJE(GeV). 

In selecting events containing r pairs with an additional radiated photon, 
background is expected from e + e~, multihadron, and two-photon events 

with any final state. Lepton pair events are selected by standard cuts |19[ against 
Z° — > qq (cut on track and cluster multiplicity), two-photon processes (cut on 
visible energy) and cosmic ray background. The cut on the acollinearity of the 
r pair is of course omitted in the preselection since it would also reject most 
of the signal events. Then e + e~ and events are recognized and rejected 

by high detected energy in the electromagnetic calorimeter (ECAL) or by high 
momentum tracks, low energy deposits in the calorimeters and /1 chamber hits, 
respectively. The observed r decay products are required to lie in a cone of half 
opening angle of 35°. We assume the r flight direction to be identical with the 
cone axis, defined by the vector sum of the associated tracks and all neutral 
clusters. A r + r~7 candidate event is selected if a photon candidate of at least 
5 GeV and less than 42 GeV energy deposit is found outside both cones. The high 
energy cut is imposed to avoid the energy region where the m T = assumption 
in the F 2 MC has an impact on the E 1 distribution. To avoid losing r decay 
products inside the beam pipe, | cos 6 T \ < 0.9 has to be valid for both r cones. 
A sample of 3 435 events survive this preselection. 
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Background from non-r events is further reduced by the following require- 
ments: 

• to suppress initial state radiation the photon has to be in the barrel region 
of the detector (| cos 6* 7 | < 0.78). 

• to reject e + e~ and fi + fi~ events, the visible energy or visible momentum of 
the more energetic r candidate must be smaller than 35 GeV. 

• the scalar sum of the momenta of the detected decay products of both 
tau candidates and the photon must be smaller than 75 GeV; furthermore 
events for which both r cones are identified as r — > [ivv decays are rejected. 
These cuts add to the suppression of events. 

• only events with 2 or 4 charged tracks (1-1 and 1-3 topologies) are retained. 

The three-body final state of the signal process e + e~ — > r + r _ 7 is completely 
determined by two independent variables, i.e. the acollinearity angle between the 
r leptons can be calculated from the measured photon energy E 1 and the mea- 
sured angle between the photon and the r~ . The measured acollinearity angle 
is required to agree within ±50° with the calculated angle. This cut greatly re- 
duces multihadron background, two-photon events and incorrectly reconstructed 
events. 

The above selection results in a total of 1429 e + e~ — > t + t~7 candidate events. 
The contribution of background from e + e~, and multihadron events to this 

sample is estimated to total (0.13 ± 0.13)%. 

Figure [| shows the observed distribution of the measured photon energy for 
the selected events. Superimposed are the expectation from the Standard Model 
(KORALZ), normalized to the number of data events (shaded) and the distribu- 
tion of F 2 -produced photons including full detector simulation (open histogram). 
In order to extract a limit on the F% form factor the data distribution of fig. |2] is 
fitted to a sum of both MC contributions using a binned likelihood £ assuming a 
Poisson distribution of the data events. In this fit the sum of both contributions 
has been normalized to the number of observed data events for each assumed 
value of F 2 . To test the method, we have performed fits to Monte Carlo event 
samples of the size of the data sample, using a 5% and a 10% F 2 contribution, 
respectively. In both cases the input values were reproduced (0.045 ± 0.010 and 
0.096 ± 0.007, respectively). 
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E y (GeV) 

Figure 2: Photon energy spectra for data (points), KORALZ (shaded histogram, 
normalized to the data), and signal (open histogram, arbitrarily normalized) . The 
dashed line shows the Monte Carlo prediction for F 2 = 0.064, also normalized to 
the data. 

Fit Results 

Figure|3|(a) shows the dependence of the likelihood function on \F 2 \. The most 
probable value is I-F2I — 0.037 which is offset from, but consistent with zero within 
about one standard deviation (^o*02s) wn i cn i s evident from the shallowness of the 
maximum of log£ in fig.||](a). The 95% confidence level value is obtained from 
fig. |3](a) at the point where log£ has dropped by 1.92 units from its maximum as 

|F 2 |< 0.064 at 95%C.L. (5) 

The analysis has also been performed by normalizing the KORALZ MC to the 
integrated luminosity of the data. In this case the luminosity has been inferred 
from the total number of r-pair events using the standard OPAL T-pair selection. 
While intuitively one might expect a tighter constraint on ji^l, the necessity to 
know the detection efficiency introduces an additional uncertainty not present 
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in the approach described above which is only sensitive to a difference in shape 
of the Ej distribution between data and MC. Both effects have been tested (see 
section on systematic errors). When normalizing to the integrated luminosity the 
uncertainty of ±6% in the knowledge of the detection efficiency yields a limit on 
\F 2 \ even slightly higher (0.065) than that obtained using the shape information 
only. 




Figure 3: (a) Likelihood as a function of the fit parameter \F 2 \. (b) Likelihood 
as a function of F 2 taking into account the interference term as described in the 
corresponding section of the text. 
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Systematic Errors 



The systematic uncertainties due to the selection cuts, the photon detection ef- 
ficiency, non-r background, binning effects, Monte Carlo statistics and normal- 
ization, and the calibration uncertainty in the photon energy measurement have 
been studied. The omission of the interference term is discussed separately below. 

Variations of the selection cuts indicate a systematic uncertainty on the limit 
on \F 2 \ of about 0.005, the largest influence coming from varying the visible 
energy cut for the r cone from 35 GeV to 32 GeV. 

The effect of binning has been studied by varying the number of bins in the 
energy spectrum from 15 to 20 and by moving the bin boundaries by half a bin 
width. This leads to a maximum change in the limit on |i*2| of +0.002. 

The photon energy calibration has been investigated as a source of systematic 
error. The agreement in the energy measurement of the electromagnetic calori- 
meter between data and MC is better than 0.9%, determined from a comparison 
of 7T° invariant masses involving all photon energies. A systematic shift of the 
photon energy by this amount results in a systematic uncertainty on the limit on 
F 2 smaller than 0.001. 

The uncertainty in the description of the photon detection efficiency and its 
dependence on the energy can make an important contribution to the systematic 
error of this analysis. An imperfect description of the efficiency would distort 
the photon energy spectrum and could thus lead to a bias for the resulting F 2 
contribution. The quality of the efficiency simulation has been checked using the 
photon energy spectrum of radiative e + e~ events where a high energy electron 
(> 43 GeV) has been required to tag the event in comparison with corresponding 
MC events. The efficiency ratio between data and MC is constant as a function of 
E 7 and consistent with unity to ±6%. The resulting effect on the limit amounts 
to less than ±0.0005. 

Background from non-r events has been estimated using the MC simulation 
considering all background reactions mentioned above and is found to be very 
small (0.13%). The total predicted background from e + e~, and hadronic 

processes amounts to 1.8 ± 1.8 events. The worst case assumption is that the 
background is distributed as the Standard Model expectation thereby artificially 
improving the limit. The resulting upper limit on F 2 , when the background is 
included, is however unchanged. 

Because the MC event sample is about 4.5 times larger than the data we have 
neglected the statistical error of the MC in the fit. To check the validity of this 
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assumption, the fit has also been performed using a method |2D which allows for 
the inclusion of both data and MC error in the likelihood. The resulting limit 
changes negligibly by 0.0003 with respect to that obtained without using the MC 
error. 

Assuming the systematic errors to be independent, they have been added in 
quadrature. Then the total systematic uncertainty has been quadratically added 
to the statistical error in each bin of the distribution and a new limit has been 
derived [] as 

|F 2 |< 0.067 at 95%C.L. (6) 



Inclusion of the Interference Term 

It has been shown [^] for the cross section of the process e + e~ — > t + t~7 that 
the contribution of the interference between the standard part and the magnetic 
part of the electromagnetic current (see eq. (|3])) can be important even if F 2 is as 
large as several percent. The authors of J21] have computed the differential cross 
section in the photon variables E 1 and cos 6* 7 for the e + e~ — > t + t~j process 
including the interference term. While this computation cannot serve for event 
simulation by means of 4- vector generation, it has been used to study the relative 
importance of both terms. 

Figure f| shows the contribution of the interference term to the cross section 
assuming ji^l = 0.064 for two different cuts on the event topology: without any 
requirement on the angular separation between the photon and the r leptons 
(thin curves), and requiring a minimum angular separation of 35° (thick curves) 
as in this analysis. In each case, the central (solid) curve shows the cross section 
assuming an I-F2I 2 contribution only, while the upper/lower (dashed) curves are 
obtained by including the interference term with a positive/negative sign. It is 
evident that the effect of the interference term is much reduced by requiring a 
minimum angular separation between the photon and the r leptons. Nevertheless, 
the effect of the interference term can be taken into account to obtain a limit on 
F 2 respecting its sign. 

A correction of the signal spectrum is obtained by reweighting the signal E 1 
spectrum according to fig. (| using different weighting factors for each value of F 2 . 
Obviously, adding the interference term with a positive sign leads to a shift of 
the photon energy spectrum towards higher energies leading to an even better 
distinction between the SM and the F 2 spectrum, while for the negative sign the 

1 Because the Poisson-based likelihood method does not have an explicit error term, x 2 has 
been used for this estimation. 
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Figure 4: Influence of the interference term; the thin lines show the cross-section 
without any additional cuts applied (solid for the \F 2 \ 2 term alone and dashed with 
interference term included). The thick lines (solid and dashed, respectively) show 
the result when a minimum angle of 35° between the photon and the r-cone axis 
is required. 

opposite is true. This treatment assumes equal efficiencies for events due to the 
interference term and to the quadratic term, an assumption which is not entirely 
correct. However, as long as the efficiency for the interference term is not larger 
than that of the quadratic term, this assumption is conservative and is, in fact, 
supported by the angular distributions shown in fig.ffl(c). The interference term 
distribution must lie in between KORALZ and the \F 2 \ 2 spectrum and due to the 
angular cuts its acceptance is smaller acceptance than that of the \F 2 \ 2 term. One 
then obtains the following 95% CL. limit using the likelihood curve of fig.^(b) 
and including systematic errors 



- 0.068 < F 2 < 0.065 



(7) 
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Discussion and Conclusions 



We have studied the reaction e + e~ — > t + t^^ to search for a contribution from 
the anomalous magnetic form factor F 2 that is related to the anomalous magnetic 
moment a T of the r lepton. The contribution of the F 2 form factor changes the 
distributions of the kinematic variables of the final state, most notably the photon 
energy spectrum. No significant contribution in addition to the Standard Model 
prediction is needed to describe the data. Comparing the data to the Standard 
Model prediction, a 95% confidence level limit on the magnitude of the magnetic 
form factor F 2 of 

\F 2 \ < 0.067 (8) 

has been placed. Taking into account the effect of the interference term between 
the Standard Model amplitude and the F 2 amplitude the 95% CL boundary on 
F 2 is 

- 0.068 < F 2 < 0.065 . (9) 



F 2 F% 

Substituting ► — the bounds on F 2 translate to limits on F 3 , the electric 

2m T e 

dipole form factor of the r lepton, for which one obtains^ 



- 3.8 x 10~ 16 e cm < eF 3 < 3.6 x 10~ 16 e cm , (10) 

with the same interpretation restrictions as mentioned for F 2 in the introduction 
and neglecting a possible influence of the r polarization on the term linear in F3. 



Appendix 



The formulae for the differential cross section |0| for e + e~ — > t + t~7 are given 
below using the Fx (SM) and the F 2 terms in the amplitude of eq. (||), but no 
interference. 



Pi 
Pi 



^ (1,0,0,1) 



4 (l,0,0,-i; 



Ptot 



pi + P2 = V« (i,o. 0,0) 



2 A compilation of recent bounds on electric and weak dipole moments of the r lepton can 
be found in 0,121. 
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q = E 7 (I,sm0 7 cos0 7 ,sin0 7 sin0 7 ,cos0 7 ) 

k\ = E T - (1, £ sin # T - 7 cos r - + ysin# T - 7 sin0 T - + icos# r - 7 ) 

&2 = Ptot -q-h 

Q M = (fci • g) fc 2 " - (fci • k 2 ) q» + (k 2 ■ q) 



x 



- cos 9 1 cos 7 

- cos Q 1 sin 7 

sin# 7 





( 


sin 7 ^ 




( 


sin 9 1 cos 7 


\ 


1- 




— COS 7 


z = q = 




sin 6* 7 sin 7 








o J 




\ 


cos# 7 


/ 



da (e + e 



ciEV- ciE^d cos 9 7 d(f> T - d(j) 7 
\ (c v 2 + c a 2 ) 



2tc 2 sin 4 9 W cos 4 9 W (s - M z o 2 ) 2 + (M z oT z u) 2 



1 [(hk 2 )(p lP2 ) - (k lPl )(k lP2 ) - (k 2Pl )(k 2P2 )} + ^ + (h<l) 



H~ 4 C^; c a 



+ F 2 2 



(kiq)(k 2 q) 

(kip 2 )(k 2 pi) - (hpi)(k 2 p 2 ) 
(hq)(k 2 q) 

a 3 



(k 2 q) (hq) 



+ q(p2- Pi) 



1 



.(hq) (k 2 q) / 



ir 2 m T 2 sin 4 9 W cos 4 8 W s 

Uc 



c') ( (p . ig) ! P2Q) + (fab)to) 



(hq)(k 2 q) 

2 



+ 2 c,, 2 (c 2 + c a 2 ) [(/cip 2 )(A;2Pi) + (fcipiXte)] 
+ 4 [(/ciP 2 )(^2Pi) - (&iPi)(&2P2)] f 
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